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. Axions are likely to be the Dark Matter (DM) that cosmological data require. They 

arise in the Peccei-Quinn solution of the strong-CP problem. In a previous work we 
showed that their model has a simple and natural generalization which yields also Dark 
Energy (DE), in fair proportions, without tuning any parameter: DM and DE arise from 
a single scalar field and are weakly coupled in the present era. In this paper we extend 
■ the analysis of this dual-axion cosmology and fit it to WMAP data, by using a Markov 

chain technique. We find that ACDM, dynamical DE with a SUGRA potential, DE with 
. a SUGRA potential and a constant DE-DM coupling, as well as the dual-axion model 

with a SUGRA potential, fit data with a similar accuracy. The best-fit parameters are 
\/-} [ however fairly different, although consistency is mostly recovered at the 2-a level. A 



peculiarity of the dual-axion model with SUGRA potential is to cause more stringent 
constraints on most parameters and to favor high values of the Hubble parameter. 
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"-H. 1. Introduction 

X 

Models with density parameters f2 ,<2e — 0.7, 0jm ~ 0.3, ~ 0.04 (for dark energy, whole 
non-relativistic matter and baryons, respectively), Hubble parameter h ~ 0.7 (in units of 100 
km/s/Mpc) and primeval spectral index n s ~ 1 fit most cosmological data, including Cosmic 
Microwave Background (CMB) anisotropies, large scale structure, as well as data on SNIa (Tegmark 
et al. 2001, De Bernardis et al. 2000, Hanany et al. 2000, Halverson et al. 2001, Spergel et al. 2003, 
Percival et al. 2002, Efstathiou et al. 2002, Riess et al. 1988, Perlmutter et al. 1988). The success 
of such ACDM models, also dubbed cosmic concordance models, does not hide their uneasiness. 
The parameters of standard CDM are still to be increased by one, in order to tune Dark Energy 
(DE). Furthermore, if DE is ascribed to vacuum, this turns out to be quite a fine tuning. 

This conceptual problem was eased by dynamical DE models (Wetterich 1988, 1995; Ratra & 
Peebles 1988, RP hereafter). They postulate the existence of an ad-hoc scalar field, self-interacting 
through a suitable effective potential, which depends on a further parameter. In RP and SUGRA 
(see below) models, this is an energy scale A (or an exponent a). 



- 2 - 



Within the frame of dynamical DE models, Mainini k, Bonometto (2004, MB hereafter) tried 
to take a step forward. Instead of invoking an ad-hoc interaction, they refer to the field introduced 
by Peccei &; Quinn (1977, PQ hereafter) to solve the strong-CP problem. If suitably tuned, such 
scheme was already shown to yield DM (Preskill, Wise & Wilczek 1983, Abbott & Sikivie 1983, 
Dine & Fischler 1983). MB slightly modify the PQ scheme, replacing the Nambu-Goldstone (NG) 
potential introduced ad-hoc, by a potential admitting a tracker solution. This scheme solves the 
strong-CP problem even more efficiently than the original PQ model. The A parameter of the 
tracker potential takes the place of the PQ energy scale, Fpq. Fixing it in the range solving the 
strong-CP problem yields DM and DE in fair proportions. Here, we shall call this cosmology dual- 
axion model. This model has several advantages both in respect to ACDM and ordinary dynamical 
DE: (i) it requires no fine tuning; (ii) it adds no parameter to the standard PQ scheme, which yields 
just DM; (iii) it introduces no field or interaction, besides those required by particle physics. This 
scheme, however, leads to predictions (slightly) different from ACDM, for a number of observables. 
In principle, therefore, it can be falsified by data. 

The essential peculiarity of the model is that it predicts a coupling between DM and DE. 
Coupled DE models were introduced by a number of authors (see, e.g., Amendola 2000, 2003; 
Gasperini, Piazza & Veneziano 2002; Perrotta & Baccigalupi 2002) as no direct evidence exists that 
DM particles follow geodesies. Most such models, however, introduce a further coupling parameter 
(3. Its tuning fixes DM-DE coupling within an acceptable range. For instance, Amendola & 
Quercellini (2003) give limits on (3 deduced from a fit to WMAP data (Spergel et al 2003); Maccio 
et al. (2004) restricted (3 even more, by studying the halo profiles, produced in N-body simulations 
with coupled DE. 

At variance from these pictures, the dual-axion scheme has no extra coupling parameter. The 
strength of the coupling is set by theory and, if this conflicts with data, the whole scheme is falsified. 
The only degree of freedom still allowed is the choice of the tracker potential. This freedom exists 
for any dynamical DE model. Also PQ exploited it by choosing a NG potential. For the sake of 
definiteness, up to now, the dual-axion scheme has been explored just in association with a SUGRA 
potential. MB showed that, in this case, the dual axion scheme predicts a fair growth of density 
fluctuations, so granting a viable picture for the Large Scale Structure. 

The recent detailed WMAP data on CMB anisotropies allow to submit the dual-axion model to 
further stringent tests, by comparing it with other cosmologies as ACDM, standard and interacting 
dynamical DE. This is done here by using a multiparameter Markov chain technique (see, e.g., 
Kosowski et al. 2003; Christensen et al. 2001; Knox et al. 2002; Lewis et al. 2002). The results 
of this numerical approach will then be discussed and understood on physical bases. None of the 
above models however performs neatly better than the others. Apparently, the best fit is obtained 
by dynamical DE based on a SUGRA potential, but its success is strictly marginal. The analysis 
of these models against CMB data is the main aim of this paper. 

The plan of the paper is as follows: In Section 2 we summarize the particle physics background 
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to the dual axion model. In particular, starting from Section 2.3, we restrict our analysis to a 
particular form of DE potential, the SUGRA potential; using a different potential might change 
some quantitative results. In Section 3 we describe the technique used to compare different models 
with CMB data and present (subsection 3.2) the results of this comparison. Section 4 is devoted 
to a final discussion of the results. 



2. A single scalar field to account for DM and DE 

Let us first remind that the strong CP problem arises from the existence in QCD (quantum- 
chromo-dynamics) of multiple vacuum states. The set of the gauge transformations £l(xi) that join 
vacuum configurations can be subdivided in classes Q n (xi), characterized by an integer n (Jackiw 
& Rebbi 1976), setting their different asymptotic behaviors. Within each class, transformations 
can be distorted into each other with continuity, while this is impossible if they pertain to different 
classes. 

Accordingly, in classical field theory there is no communication between different-n gauge 
sectors. In quantum field theory, instead, tunneling is possible thanks to instanton effects, so 
that any vacuum state is a superposition of the vacua |0 n ) (of the nth sector), of the kind \0g) = 
En |0„)exp(m6»). 

The effects of varying the ^-vacuum can be recast into variations of a non-perturbative term 

£ e = ^-9G-G (1) 

(a s : strong coupling constant, G and G: gluon field tensor and its dual) in the QCD Lagrangian 
density. However, chiral transformations also change the vacuum angle, so that the ^-parameter 
receives another contribution, arising from the EW (electro-weak) sector, when the quark mass 
matrix M. is diagonalized, becoming 

e eff = 6 + Arg det M . (2) 

The Lagrangian term (1) can be reset in the form of a 4-divergence and causes no change of the 
equations of motion. It however violates CP and, among various effects, yields a neutron electric 
moment d n ~ 5 ■ 10~ 16 # e jj e cm, conflicting with the experimental limit d n < 10~ 25 e cm, unless 
*e//<10- 10 . 

The point is that the two contributions to 6 e ff are uncorrelated, so that there is no reason 
why their sum should be so small. 

PQ succeed in suppressing this term by imposing an additional global chiral symmetry U(1)pq, 
spontaneously broken at a suitable scale Fpq. The axion field is suitably coupled to the quark 
sector. The details of this coupling depend on the model and may require the introduction of an 
ad-hoc heavy quark (Kim 1979, Shifman et al. 1979, see also Dine et al. 1981, Zhitnisky 1980). 
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The U(1)pq symmetry suffers from a chiral anomaly, so the axion acquires a tiny mass because of 
non-perturbative effects, whose size has a rapid increase around the quark-hadron transition scale 
h-QCD- The anomaly manifests itself when a chiral U(1)pq transformation is performed on the 
axion field, giving rise to a term of the same form of (1), which provides a potential for the axion 
field. 

As a result, 9 is effectively replaced by the dynamical axion field. Its oscillations about the 
potential minimum yield axions. This mechanism works independently of the scale Fpq. Limits 
on it arise from astrophysics and cosmology, requiring that \Q l0 GeV < Fpq < 10 12 GeV; in turn, 
this yields an axion mass which lays today in the interval W~ 6 eV < rriA < 10~ 3 eV. 

More in detail, in most axion models, the PQ symmetry breaking occurs when a complex scalar 
field $ = (f>e id /V2, falling into one of the degenerated minima of a NG potential 

V(<S>) = \M 2 -F 2 PQ f , (3) 

develops a vacuum expectation value {(f)) = Fpq. 

The CP- violating term, arising around quark-hadron transition when qq condensates break 
the chiral symmetry, reads 



V, = 



^2(Q(T)\qq\0(T))m q 



(l-cos0) (4) 



(So extends over all quarks), so that 9 is no longer arbitrary, but shall be ruled by a suitable 
equation of motion. The term in square brackets, at T ~ 0, approaches m 2 f 2 (m, and f n : ir- 
meson mass and decay constant). In this limit, for 9 <C 1 and using A = 9 Fpq as axion field, 
eq. (4) reads: 



V 1 ~ tq\m q )mlfl^- ■ (5) 



1 2 / \ 2 c2 

-q KJmJ,^- 
1 *PQ 

here q{m q ) is a function of the quark masses qf, in the limit of 2 light quarks (u and d), q = 
y/ m u /md{^--\-m u /m l i)~ l ■ The A field bears the right dimensions but, here below, will no longer be 
used, and the axion degrees of freedom will be described through 9 itself. Eq. (5), however, shows 
that, when (qq) is no longer zero (since T < Aqcd), the axion mass decreases with temperature 
approaching the constant value ijia = m n f w q(m q )/FpQ for T <C Aqcd- 

Accordingly, the equation of motion, in the small 9 limit, reads 

9 + 2-9 + a 2 m\9 = , (6) 
a 

(here a is the scale factor and dots yield differentiation with respect to conformal time, see next 
Section), so that the axion field undergoes (nearly) harmonic oscillations, as soon as ijia exceeds 
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the expansion rate; then, his mean pressure vanishes (Dine & Fischler 1983) leaving axion as a 
viable candidate for cold DM. 

MB replace the NG potential in eq. (3) by a potential V($) admitting a tracker solution 
(Wetterich 1988,1995, RP 1988, Ferreira & Joyce 1998, Brax & Martin 1999, 2001, Brax, Martin 
& Riazuelo 2000). The field $ is complex and V($) is £7(1) invariant, but there is no transition to 
a constant value Fpq, which is replaced by the modulus 4> itself, slowly evolving over cosmological 
times. At a suitable early time, quantum dynamics starts to be fairly accounted by the potential 
V; soon after, 4> settles on the tracker solution in almost any horizon, however breaking the U(l) 
symmetry, by the values assumed by 8, in different horizons. Later on, when chiral symmetry 
breaks, dynamics becomes relevant also for the 8 degree of freedom, as in the PQ case. At variance 
from it, however, this happens while <p continues its slow evolution, down to the present epoch, when 
it accounts for DE. Owing to the <j> evolution, however, the axion mass evolves, over cosmological 
times, also for T <C h-QCD (see below). 

The <3? field, therefore, besides of providing DM through its phase 9, whose dynamics solves 
the strong CP problem, also accounts for DE through its modulus (f>. 

This scheme holds for any DE potential admitting tracker solutions. To be more specific, MB 
use a SUGRA potential (Brax, Martin k Riazuelo 2000; Brax & Martin 1999; Brax & Martin 2001) 
finding that, at the quark-hadron transition, <f> can be naturally led to have values ~ Fpq, increasing 
up to ~ m p = G~ x l 2 (the Planck mass), when approaching today. The only free parameter is the 
energy scale in the SUGRA potential, that must be ~ 10 10 GeV. With this choice, 9 is driven to 
values even smaller than in the PQ case, so that CP is apparently conserved in strong interactions, 
while f2 ,m, Q ,de (the DE density parameter) and £l ^ take fair values. 



2.1. Lagrangian theory 

In the dual-axion model we start from the Lagrangian 

£ = V^g{g^d^d^ - v(<f>)} , (7) 

which can be rewritten in terms of <p and 9, adding also the term breaking the U(l) symmetry. 
Then it reads: 

C = ^\g^4>du^ + ^0^98,9] -V(^) - m 2 (T,<P)^\l - cos9)^ . (8) 

Here g^ u is the metric tensor. We shall assume that ds 2 = g^dx^dx" = a 2 (dr 2 — rjijdxidxj), so 
that a is the scale factor, r is the conformal time; Greek (Latin) indexes run from to 3 (1 to 3); 
dots indicate differentiation in respect to r. The mass behavior for T ~ Aqcd will be detailed in 
Section 2.2. The equations of motion, for the <j> and 9 degrees of freedom, read 

9 + 2 (^ + 1) ^ + ™Vsin# = , (9) 
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4> + 2-<j) + a 2 V'(<j)) = 4>9 2 . (10) 
a 

(Notice that m 2 (T,(j))(f) 2 is 4> independent, see below.) In what follows, the former equation will 
be always considered when sin ~ 0. In particular, taking into account the condition 0^1, the 
expressions for the energy densities p e ^ = pe, < p-kin + Pe,4>; P ot and the pressures p e ^ = pe^-kin- Pe,<j>; P ot 
are obtainable by combining the terms 



^ 2 m 2 (T, 

-if on . . i /// i .? (.}\:) i i — ens f / i ' 

2a 2 



P9,kin = —»- , Pe, P ot = m 2 (T, 0)^(1 - cos 0) ~ " ' ' 2 " 2 



</> 2 

P4>,kin = ^2 ' P4>,Pot = V (<t>) ■ ( n ) 

When underg oes many (nearly) harmonic oscillations within a Hubble time, (pofein) — {p9, P ot) 
and (jpq) vanishes (Dine & Fischler 1983). Under such condition, using eqs. (9) and (10), it is easy 
to see that . 

CL Tfl CL TYl 

pe + 3-p e = —pe , p<p + 3-(p</, + p^) = pe ■ (12) 

a m a m 

When m is given by eq. (15) and (16) here below, rn/m = —<j)/4> — 3.8T/T. At T ~ 0, instead, it 

is just rh/m ~ —(f>/<f>- 

The and <f> components account for DM and DE, respectively. Accordingly, in the sequel, 
the indexes q, <p will be replaced by dm-, de- 

Eqs. (12) clearly show that an exchange of energy occurs between DM and DE. From this point 
of view, the MB model belongs to the set of coupled models treated by Amendola (2000, 2003). It 
is however characterized by a time-dependent coupling. In fact, in the small limit, after averaging 
over cosmological times, the r.h.s. of eqs. (10) and (12) read C((p)(pQ)a 2 and ±C(0) 4> {pe)-> if we set 
C{<p) = Here C is the DE-DM coupling introduced by Amendola (2000, 2003), who however 
studies extensively only the case C = ^(Wir/Sm 2 ) 1 / 2 , with constant j3. In these latter DE models, 
a <p-MDE phase takes place after matter-radiation equivalence. It differs from a matter dominated 
expansion because of the contribution of the kinetic part of the DE field to the expansion source. 
A regime of this kind is present also in the dual-axion model and is shown in Fig. 5, here below. 
Because of the <f> dependence, however, the DM-DE coupling, in the dual-axion model, weakens as 
we approach the present cosmological epoch. 

The most stringent limits on (3 are set by non-linear predictions (Maccio et al. 2004) and 
restrict (3 to values < 0.1-0.2, in order to avoid a too high concentration in DM halos. In turn, 
this is due to the behavior of the effective mass of DM particles, in the presence of the /3-coupling. 
A preliminary inspection indicates that the behavior expected here is opposite and that halo con- 
centrations, in average, should be smaller than in ACDM. This point, however, must be inspected 
in much more detail. 

Let us also notice that the former eq. (12) can be integrated soon, yielding p& m oc m/a 3 . In 
particular, this law holds at T <C A.qcd- Accordingly, at late times 

Pdma 3 (f> - const., (13) 
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so that the usual behavior pd m oc a~ 3 is modified by the energy outflow from DM to DE. This 
modification is stronger when <p varies rapidly and is damped when <f> attains a nearly constant 
behavior. 



2.2. Axion mass 

According to eq. (9), the axion field begins to oscillate when: 

m(T, 0)a ~ 2 ( - + ^ ) . (14) 




In the dual-axion model, just as for PQ, the axion mass rapidly increases when the chiral symmetry 
is broken by the formation of the qq condensate at T ~ Aqcd- In the dual-axion model, however, 
the axion mass varies also later on, because of the evolution of </>, when m(T, <p) is 

m {<t>) = ^ ; (15) 

with ij,qcd — 80MeV. Since (p ~ m p today, the present axion mass m a ~ 5 • 10~ 13 eV. At high 
temperature, according to Gross et al. (1981), 

m(T^)^O.lm^)^) 3 ' 8 (16) 

This expression must be interpolated with eq. (15), to study the fluctuation onset for T ~ Aqcd- 
We report the results of the solution of the equations in Sec. 2.1, obtained by assuming 

(0 1 1/3 ' 8 A \ 3 ' 8< - 1_a ^ ac - ) 
T QC ° I a<a c (17) 

m(T, <f>) = m ((f>) a > a c (18) 

with a c = T / Aqcd 0.1 1 / 3 - 8 = 2.16 • 10" 12 ; here T Q = 2.35 • 10~ 4 eV, Aqcd = 200MeV~ 2.5/i Q cD- 
The expressions (17), (18) assure that m(T, (f>) and its first derivative are continuous and that the 
axion mass meets its low-T behavior m (T) at T ~ 0.55 Aqcd. The effects of selecting a different 
value for Aqcd were studied. Also the impact of the assumed T dependence was considered. A 
change of Aqcd, by a factor 2, yields a variation of A in the SUGRA potential (see next section) 
by not more than a few percent. Taking a power v ^ 1 of (1 — a/a c ), in the exponent, instead, 
does not affect continuity but changes the rapidity by which the high-T regime (16) is met. In 
Fig. 4 we plot both the interpolated mass behavior and the first 6 oscillations. Changing v just 
slightly displaces the interpolating curve between ~ oqcd and ~ 1.5 clqcd, causing a minor phase 
shift. Let us finally outline that eqs. (15) and (16), as well as the interpolation (17), (18), assure 
that m(T, <j))(f) is ^-independent, as is required to give the equation of motion the form (10). 
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2.3. Using the SUGRA potential 

Most above features are general and do not depend on the choice of the potential V. We shall 
now assume that 

(SUGRA potential; see Brax & Martin 1999, 2001, Brax, Martin k Riazuelo 2000). Apparently 
V(&) depends on the two parameters A and a. When they are independently assigned, is also 
fixed. Here we prefer to use Qde and A as independent parameters. The latter scale is related to a 
as shown in Figure 1. This relation is almost independent from the presence of DE-DM coupling. 
The slight shifts due to the couplings, evaluated through the expressions of this paper, are just 
slightly above the numerical noise and are also shown in Figure 1. 

The potential (19) does not depend on 9 and, in the radiation dominated era, admits the 
tracker solution 

0«+ 2 = g a A a+4 a 2 T 2 , (20) 

with g a = a(a + 2) 2 /4(a + 6). This tracker solution, characterizing SUGRA models at very high 
z, is abandoned, because of the DE-DM coupling, when the term <f>9 exceeds a 2 V, and the field 
enters a different tracking regime: 

<P 2 = \p dm a 2 T 2 . (21) 

Fig. 2 is a landscape behavior of densities, starting from the high-z tracking regime (20), passing 
then to the new intermediate tracking regime, and reaching the regimes when DE density eventually 
exceeds first radiation (z ~ 100), then baryons (at z ~ 10) and DM (at z ~ 3). In Fig. 3, we show 
the evolution of the DE-field 4>. In Fig. 4, we focus on the Q-H transition and magnify the scale 
dependence of the mass and the onset of 9 oscillations. Finally, in Fig. 5, we show the behaviors of 
the density parameters fij (i = r, b, 9, (/), i.e. radiation, baryons, DM, DE) vs. the scale factor a. 
More detailed pictures of the single transitions are shown in MB. 



2.4. Parameter fixing 

In general, once f2 <fe and h are given, a model with dynamical (coupled or uncoupled) DE is 
not yet determined, as a or A are still to be fixed. When fitting WMAP data, here below, we shall 
mostly refer to the energy scale A and remind that, when £l 0) de and A are given also a is fixed. 

Other potentials show similar features and this implies that, in general, dynamical DE models 
depend on one more parameter than, e.g., ACDM; accordingly, one can expect that they more 
easily accommodate observational data. In principle, data fitting is even more facilitated by DM- 
DE coupling, which adds a further parameter (the strength of coupling) to the theoretical plot. 

In the dual-axion model with a SUGRA potential, no such arbitrariness exists. If this model, 
where a suitable DM-DE coupling is present, succeeds in fitting observational data, this will not be 
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favored by extra parameters. The coupling, in particular, depends on no parameter and the very 
scale A is set by theoretical consistency arguments, similarly to what happens for the Fpq scale in 
the PQ approach. 

Let us follow the behavior of pdm, backwards in time, until the very beginning of the oscillatory 
regime, when the approximation S < 1 begins to hold. Using the scaling law (13) until then and, 
at earlier times, the more general laws (11), together with eq. (14), we build a system yielding the 
scale factor when the fluctuations start, and the scale A in the SUGRA potential (19). 

It turns out that ah and A are almost independent from £l ,dm and h, and we shall now outline 
the qualitative reasons of this feature. According to eqs. (11) and (15), at the ah scale it must be 

/ \ 2 

Ph,dm - m 2 h <p 2 {a h ) e 2 (a h ) = ( — J p A QCD e 2 {a h ) , (22) 

where rrih = m(a^), while e 2 = 2((1 — cos 6)) is ~ 6 2 when the oscillatory regime is onset. Owing to 
eq. (14), however, the oscillation onset occurs when the axion mass is approximately the inverse of 
a h T ( a h) an d the latter quantity can be related to the temperature Th, via the Friedmann equation. 
This yields that T% ~ m p mh/8 and, using again eq. (15), that 

T h ^5.10' 2 p QCD (^^-) 1/2 • (23) 



h m J 

Let us now assume that it is now (f> ~ m p and use the law (13) from now to ah- It is then 



Th ~ hqcd 



5 • 10 3 



provided that fi QC D ^ 4 • 10 n T o and p o4m ~ 3 • 10 4 T o 4 O Ojdm /i 2 . 



(24) 



The temperature Th (and the scale ah) are then essentially model independent because of the 
power 1/5 at the l.h.s. of eq. (24). If we consider the plots in Fig. 4, obtained through a numerical 
integration for fl 0) d m = 0.27 and h = 0.7, we see that (m^/m ) 3 e 2 ~ 10~ 2 -10~ 3 , so that Th falls 
around Aqcd- However, no appreciable displacement can be expected just varying Qo,dm and 
h, even though we suppose that this induces substantial variations of the factor (mh/m ) 3 e 2 (a^) 
which, however, cannot exceed unity or lay below ~ 1CP 3 . 

The model independence of ah implies that the scale A is also almost model independent. The 
numerical result ah — 1CP 13 is consistent with A ~ 1.5 • 10 10 GeV, but neither value will change 
much just by varying 0, o ,dm- In practice, when £l ,dm goes from 0.2 to 0.4, log 10 (A/GeV) (almost) 
linearly runs from 10.05 to 10.39 and ah steadily lays at the eve of the quark-hadron transition. 

The residual dependence of a and A on fi 0)( / m is plotted in Fig. 6, together with the correspond- 
ing values for w at z = 0. Significantly smaller values for A are obtainable only for unphysically 
small DM densities. 
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The only way to modify the result is to consider values of <p significantly different from m p , 
today. In fact, the rest of the plot can be understood as an effect of the position reached by the 
(j) field, in the SUGRA potential, at the present time. The peak of w corresponds to a maximum 
of kinetic energy and occurs if the minimum of SUGRA potential is attained today. Today lays in 
the proximity of the minimum for a fairly wide A interval. Accordingly, a significant A variation, 
thereabout, yields just a modest shift of Q ,dm, as is shown in the lower plot. For still greater 
A's, (j) would still be in its pre-minimum descent. For even smaller amounts of DE, the present (j) 
configuration would still lay in the <j)-MD era, when kinetic energy dominates for DE, although the 
very DE contribution to the overall energy density becomes negligible. 

A model with DE and DM given by a single complex field, based on a SUGRA potential, 
therefore bears a precise prediction on the scale A, for the observational £l ,dm range. Moreover, 
just the observational £l ,dm range (0.2-0.4) corresponds to <fi ~ m p today. 

The rest of this paper is devoted to a comparison of this and other models against WMAP 
data. In such comparison, however, the scale A shall not be fixed, a priori, in the expected range. 
In general, the dual axion model belongs to a wider class of coupled DM-DE models, which can be 
dubbed ^ _1 -models because of the shape of the coupling, where the scale A is left free. 

The features of 4>~ l models, in general, are quite reasonable. They have two independent dark 
components, DM and DE, suitably coupled, about whose nature no assumption is made, just as 
in standard coupled DE models. In the latter models, however, the strength of the coupling is 
gauged by an extra parameter (5. <f)~ l models, from this point of view, already offer less freedom, 
as C = <p~ l and no coupling modulation is allowed. 

Should the comparison of (j)" 1 models with WMAP data favor parameter values compatible 
with cosmological observables and values of A compatible with the dual-axion model, this means 
that WMAP data support the origin of DE and DM that the dual-axion model suggests. 

3. Comparison with WMAP data 

WMAP data have been extensively used to provide tight constraints on cosmological parame- 
ters. They consists of high precision estimates of the anisotropy power spectrum Cj up to I ~ 900, 
as well of the TE correlation power spectrum Cf E up to I ~ 450. We shall use these data to con- 
strain possible cosmologies, in a parameter space of 7 to 8 dimensions. A grid-based likelihood anal- 
ysis would then require prohibitive amounts of CPU time and we use a Markov Chain Monte Carlo 
(MCMC) approach, as it has become customary for CMB analysis (e.g., Christensen et al. 2001; 
Knox et al. 2002; Lewis et al. 2002; Kosowski et al. 2002; Dunkley et al. 2004). 

The principal analysis of WMAP first-year data (Spergel et al. 2003) constrained flat ACDM 
models defined by six parameters: Vt ^h 2 , fJ 0j?n /j 2 , h, n, the fluctuation amplitude A and the optical 
depth r. Notice that, with the naming convention used here, Ojm = f2 0i 6 + Q ,dm- As possible 
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extensions of ACDM cosmologies, several works considered models with a fixed state parameter 
w = Pde/Pde (e.g., Spergel et al. 2003; Bean & Dore 2004; Tegmark et al. 2004; Melchiorri 2004), 
or adopted z-dependent parameterizations of w{z) interpolating between early-time and late-time 
values (e.g., Corasaniti et al. 2004; Jassal et al. 2004; Rapetti et al. 2004). A general conclusion 
was that current data mostly allow to constrain only the present state parameter, w(z = 0) < —0.80. 

In this work we consider, instead, three classes of dynamical DE, requiring the introduction of 
additional parameters specifying the physical properties of the scalar field, (i) SUGRA dynamical 
DE require the introduction of A = log 10 (A/GeV), yielding the energy scale in the potential (12). 
(ii) In constant coupling DE, the coupling parameter (3 = C (Srn^/lQit) 1 / 2 is also needed, (iii) In 
the case of the dual-axion model, the last parameter is excluded, for it is simply C = 4>~ l . Also 
the scale A is constrained by the requirement that f2 ,de lays m a feh range (also solving the strong 
CP problem). Hence, in the dual-axion model, A and £l ,de axe no longer independent parameters. 
We however consider a wider class of coupled DE models, that we call models, leaving A as a 
free parameter. Our aim is to test whether WMAP data constrain it into the region turning a 
model into a dual-axion model. 

In the use of MCMC, as well as in any attempt to fit CMB data to models, a linear code provid- 
ing Q's is needed. Here we use our optimized extension of CMBFAST (Seljak & Zaldarriaga 1996), 
able to inspect the cosmologies (i), (ii) and (iii). Then, the likelihood of each model is evaluated 
through the publicly available code by the WMAP team (Verde et al. 2003) and accompanying 
data (Hinshaw et al. 2003; Kogut et al. 2003). 

3.1. Implementing a MCMC algorithm 

A MCMC algorithm samples a known distribution £(x) by means of an arbitrary trial distri- 
bution p(x). Here £ is a likelihood and x is a point in the parameter space. The chain is started 
from a random position x and moves to a new position x', according to the trial distribution. The 
probability of accepting the new point is given by £(x')/£(x); if the new point is accepted, it is 
added to the chain and used as the starting position for a new step. If x' is rejected, a replica of x 
is added to the chain and a new x' is tested. 

In the limit of infinitely long chains, the distribution of points sampled by a MCMC describes 
the underlying statistical process. Real chains, however, are finite and convergence criteria are 
critical. Moreover, a chain must be required to fully explore the high probability region in the 
parameter space. Statistical properties estimated using a chain which has yet to achieve good con- 
vergence or mixing may be misleading. Several methods exist to diagnose mixing and convergence, 
involving either single long chains or multiple chains starting from well separated points in the 
parameter space, as the one used here. Once a chain passes convergence tests, it is an accurate 
representation of the underlying distribution. 

In order to ensure mixing, we run six chains of ~ 30000 points each, for each model cate- 
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gory. We diagnose convergence by requiring that, for each parameter, the variances both of the 
single chains and of the whole set of chains (W and B, respectively) satisfy the Gelman & Rubin 
test (Verde et al. 2003; Gelman & Rubin 1992), R < 1.1 with: 

R = [{N-l)/N]W + {l + l/M)B _ (25) 

Here each chain has 2N points, but only the last N points are used to estimate variances, and 
M is the total number of chains. In most model categories considered, we find that the slowest 
parameter to converge is A. 



3.2. Results 

The results of this paper mostly concern the fit of the double-axion model with WMAP data. 
Let us however also briefly outline the features of the model. 

Quite in general, cosmological models fitting observations require a a triple coincidence between 
DM, baryon and DE densities, the last one occurring just at the present epoch. Any such coincidence 
requires the tuning of a suitable parameter. It is then natural to look for an underlying physics, 
able to predict the parameter values that cosmology requires. 

If DM is made of axions, we must tune the Fpq parameter in the NG potential involving a $ 
field, whose phase 6 is then the axion field. Here we saw that the PQ model does not strictly require 
a NG potential. A potential causing no immediate settling at an energy minimum, but inducing 
a slow rolling along a tracking solution, achieves the same aims. If such potential is SUGRA, it 
contains an energy scale A to be tuned, in stead of Fpq. If this is done, so to yield the required 
amount of DM, this scheme provides, as an extra bonus, a fair amount of DE, which is prescribed 
to interact with DM in a peculiar way, never studied in the literature. 

The possible shapes of DE-DM interaction considered up to now involved an ad-hoc coupling 
parameter. No such parameter exists here. For instance, a standard model of coupled SUGRA 
involves a density parameter £l 0i dmi an energy scale A and a coupling parameter (3. Here, not 
only f3 no longer exists, but also A is fixed once Qo : dm is assigned. The model is therefore highly 
constrained and data can easily falsify it. The only available degree of freedom, to modify it, 
amounts to replacing the SUGRA potential with another potential shape. Another reserve to be 
born in mind is that axion models have a contribution to DM coming from the decay of topological 
structures, which were not considered here. 

These arguments outline the significance of the fit of the model with WMAP data. The basic 
results of this work concern this fit and are summarized in the Tables 1-3. For each model category 
we list the expectation values of each parameter and the associated variance; we also list the values 
of the parameters of the best fitting models. The corresponding marginalized distributions are 
plotted in figures 7-9, while joint 2D confidence regions are shown in figures 10-12. 
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The values of x 2 > for each category of models can be compared, taking into account the number 
of degrees of freedom. This comparison is shown in Table 4. The smallest x 2 is obtained for the 
uncoupled SUGRA model, which performs slightly better than ACDM. Differences however are 
really small and yield no support to any model category. 

It must be however reminded that the ^> _1 models, whose fitting results are reported in Table 3 
and Figures 9 and 12, include the dual-axion model, but many other cases as well. Our approach 
was meant to test whether CMB data carry information on A and how this information fits the A 
range turning a model into the dual-axion model. 

Let us also outline that Tables 1 and 2 and the corresponding figures, concerning uncoupled 
or constant- coupling SUGRA models, outline that WMAP data provide no real constraint on 
A = log(A/GeV), when allowed to vary from ~ —12 to 16. No limitation exists even on its sign. 
On the contrary, when a (f>~ 1 coupling is set, loose but precise limitations on A arise, as is shown 
in Table 3. In the presence of this coupling, the 2-a A-interval ranges from ~ 10 to ~ 3 • 10 10 GeV, 
including the range required by the dual-axion model. 

In fact, in this case, the location of peaks in spectra is more strictly related to the 4> evolution. 
In turn, only a restricted sect of A values allows an /-dependence of multipoles consistent with 
data. Fig. 13 shows this fact for anisotropy, while Figs. 14 and 15 are predictions for the TE and 
E-polarization spectra. 

Fig. 13 also shows why no model category neatly prevails. At large I all best-fit models yield 
similar behaviors. In turn this shows that discrimination could be achieved by improving large 
angular scale observation, especially for polarization, so to reduce errors on small-/ harmonics. 

4. Discussion 

Our analysis does not concern the double axion model only, but the whole spectrum of models 
built using a SUGRA potential. This is useful also to evaluate the significance of the fit of the 
dual-axion model. 

4.1. Uncoupled and coupled SUGRA models 

A first point we must therefore outline is that SUGRA uncoupled models are consistent with 
WMAP data. The ratio w = p/p, for most these models is < —0.80 at z = 0. However, they 
exhibit a fast variation of w, which already attains values ~ —0.6 at z ~ 1-2. This sharp decrease 
does not conflict with data and these models perform even better than ACDM. For uncoupled 
or constant-coupling SUGRA models, the analysis in the presence of priors leads to analogous 
conclusions. 
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Table 1: SUGRA parameters in the absence of DE-DM coupling: for each parameter x, the expec- 
tation value (x) , variance a x , and maximum likelihood values x ma x , in the 7-dimensional parameter 
space, are shown. 
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Table 2: SUGRA parameters in the presence of a constant DE-DM coupling [3: The parameter 
space is 7-dimensional and parameter values are shown as in the previous Table. 
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Table 3: SUGRA parameters for a (ft model. The parameter A is left arbitrary; at variance 
from other model categories, A here is constrained and consistency with the dual-axion model is 
recovered at 2cr's. Parameter values are shown as in Table 1. 
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Fig. 1. — Values of a corresponding to given A scales in SUGRA models with Q ,dm = 0.27. The 
tiny differences, just above noise level, for coupled models (ap and oti/$ for constant coupling and 
(f)" 1 coupling, respectively), are also shown. 



Table 4: Goodness of fit. For each class of DE considered, Table lists the number of degrees of 
freedom (d.o.f.), the reduced Xeffi an d the corresponding probability of the best-fit model. Figures 
for ACDM and dual-axion models are also included. ACDM models have 1342 degrees of freedom, 
uncoupled and (ft -1 have 1341, while fixed-/? has 1340. 

xlff Prob. 

no coupling 1.064 5.0 % 

/3-coupling 1.066 4.7 % 

^-coupling 1.074 2.9 % 

dual-axion 1.081 2.0 % 

ACDM 1.066 4.7 % 



-16- 




-15 -10 -5 

log(a/a ) 



Fig. 2. — Behaviors of densities from the tracking regimes down to the regimes when DE density 
exceeds radiation (z ~ 100), baryons (at z ~ 10) and DM (at z ~ 3). 




log(a/a G ) 



Fig. 3. — Evolution of the DE-field. Notice the jump around the QH scale and the low-z rebounce. 
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Fig. 4. — Interpolated mass behavior and first 6 oscillations, clqcd is the scale when T = Aqqd- 
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Fig. 5. — Density parameters Qi (i = r, b, 6, 4>, i.e. radiation, baryons, DM, DE). 
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Fig. 6. — The lower plot shows the variation of the exponent a and the A scale in the SUGRA 
potential, when varying the present fraction of DM. The faster A increase, at £l ,dm — 0-6, is due 
to a stationary behavior in the set of models where the 4> fields attains (approximately) today the 
minimum of the SUGRA potential. In fact, if it were Q ,dm — 0.6 today, the upper plot shows 
that w would corresponds to top kinetic energy and minimum potential energy. These plots are 
obtained for h = 0.7. 
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Fig. 7. — Marginalized distributions for the 7-parameters SUGRA model with no priors (solid 
lines), BBNS prior (long dashed) or HST prior (dot-dashed). Short dashed (dotted) vertical lines 
show the boundaries of 68.3 % c.l. (95.4 % c.l.) interval; for A only upper limits are shown. 
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Fig. 8. — As Fig. 7 but for the 8-parameters constant coupling model. For A and (3 only the upper 
c.l. boundaries are shown. 
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Fig. 9. — As Fig. 7 but for the 7-parameters <j) 1 model. 
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Fig. 10. — Joint 2D constraints for SUGRA models. Light (dark) shaded areas delimit the region 
enclosing 68.3 % (95.4 %) of the total points. 
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Fig. 11. — As Fig. 10 but for constant coupling models. 
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Fig. 12. — As Fig. 10 but for (j) 1 models. Here, cosmological parameters are more stringently 
constrained than in other models. 
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Fig. 13. — Cf spectra for the best fit SUGRA (solid line), constant coupling (dotted line), 
coupling (dashed) and dual-axion (dot-dashed) models. Dual-axion model results from considering 
only those (f)~ l with 9.5 < A < 10.5. The binned first-year WMAP data are also plotted. 
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Fig. 14.— Best fit Cf h spectra. 
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Fig. 15. — Best fit Cf spectra. 
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Best fit cosmological parameters exhibit some dependence on the model. First, the opacity r 
is pushed to values exceeding the ACDM estimates (see also Corasaniti et al. 2004). This can be 
understood in two complementary ways: (i) Dynamical DE models, in general, exhibit a stronger 
ISW effect, as the field <f> itself varies during the expansion, and DE effects extend to greater z. This 
increases Cf in the low-/ plateau (e.g., Weller & Lewis 2003). To compensate this effect the fit 
tends to shift the primeval spectral index n s to a greater value. Owing to the r-n degeneracy, this 
is then compensated by increasing r. (ii) In dynamical DE models, the expected TE correlation, at 
low I, is smaller than in ACDM (Colombo et al. 2003). A given observed correlation level, therefore, 
requires a greater r. In any case, values of r ~ 0.07 keep consistent with data within less than 
2-cr's. 

Greater r's push upwards Qf,h 2 estimates, although best fit values are consistent with ACDM 
within 1— a. Adding a prior on Q b h 2 = 0.0214 ± 0.0020 (BBNS estimates, see, e.g. Kirkman et al. 
2003) lowers h, within 1-a from HST findings. We therefore also consider the effect of a prior on 
h. In Figures 7 and 8 the effects of priors are shown by the dashed red line (prior on f^/i 2 ) and 
the dot-dashed blue line (prior on h). 

The former prior affects mainly reionization and n s ; r and n s are lowered to match WMAP's 
findings, the high tail of r distribution is partly suppressed. The physical analysis of primeval 
objects causing reionization (e.g., Ciardi et al. 2003; Ricotti & Ostriker 2004) could however hardly 
account for values of r ~ 0.3 which are still allowed but certainly not required. 

The latter prior favors greater h values. In the absence of coupling, this favors low-A models, 
closer to ACDM. In fact, the sound horizon at decoupling is unaffected by the scale A, while the 
comoving distance to last scattering band is smaller for greater A's. Then, as A increases, lower 
h values are favored to match the angular position of the first peak. In the presence of coupling, 
there is a simultaneous effect on (3, as greater /3's yield a smaller sound horizon at recombination, 
so that the distribution on h is smoother. 

A previous analysis of WMAP limits on constant coupling models had been carried on by 
Amendola & Quercellini (2003). Their analysis concerned potentials V fulfilling the relation 
dV/d(j) = BV N , with suitable B and N. Furthermore, they assume that r = 0.17. Our analy- 
sis deals with a different potential and allows more general parameter variations. The constraints 
on P we find are less severe. It must be however outlined that (5 <> 0.1-0.2 seem however forbidden 
by a non-linear analysis of structure formation (Maccio et al. 2004). 

4.2. Dual axion models 

Let us consider then </> _1 — models, which generalize dual-axion models to an arbitrary A scale, 
used as a parameter to fit data. 

Parameters are better constrained in this case, although the overall model likelihood is similar. 
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This is made evident by Fig. 12. In particular, at variance from the former case, the energy scale A 
is significantly constrained and, within 2-cr, constraints are consistent with the double-axion model. 

Several other parameters are constrained, similarly to coupled or uncoupled SUGRA models. 
What is peculiar of (/) _1 -models is the range of favored h values: the best-fit 2-cr interval does not 
extend much below 0.85 . 

This problem is slightly more severe for the dual-axion model. The point is that this model 
naturally tends to displace the first Cf peak to greater I (smaller angular scales) as coupling does, 
in any case. But, in the absence of a specific coupling parameter, the very intensity of coupling, 
in these models, depends just on the scale A. Increasing A requires a more effective compensation, 
favoring greater values of h. 

Apart of the possibility that h is currently underestimated, let us remind other substantial 
options which are still to be deepened: (i) The contribution of topological singularities to DM were 
not considered here and taking them into account could naturally increase the amount of DM, 
yielding smaller A values for the double-axion model; fig. 6 however shows that this is hardly 
an efficient solution, (ii) The choice of a SUGRA potential is however arbitrary, the dual axion 
model does not require SUGRA. Other potentials can possibly yield the same coupling intensity in 
agreement with a smaller h. 

5. Conclusions 

The first evidences of DM date some 70 years ago, although only in the late Seventies limits 
on CMB anisotropics made evident that a non-baryonic component had to be dominant. DE could 
also be dated back to Einstein's cosmological constant, although only SNIa data revived it, soon 
followed by data on CMB and deep galaxy samples. 

Axions have been a good candidate for DM since early Eighties, although various studies, as 
well as the occurrence of the SN 1987a, strongly constrained the PQ scale around values 10 10 < 
FpQ % 10 12 GeV. Contributions to DM from topological singularities (cosmic string and walls) 
narrowed the constraints to Fpq. Full agreement on the relevance of these contributions has not 
yet been attained and, in this paper, they are still disregarded. 

The fact that DM and DE can both arise from scalar fields, just by changing the power of 
the field in effective potentials, already stimulated the work of various authors. A potential like 
(19) was considered in the so-called spintessence model (Boyle et al. 2002; Gu Sz Hwang 2001). 
According to the choice of parameters, $ was shown to behave either as DM or as DE. In the frame 
of a model of tachyon DE, Padmanabhan &; Choudhury (2002) also built a model where DM and 
DE arise from a single field. The possibility that both DM and DE arise from the solution of the 
strong CP problem was also suggested by Barr & Seckel (2001). Their model, however, does not 
deal with dynamical DE and aims to explain why the vacuum energy is so finely tuned, while DM 
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is simultaneously provided. 

Here we deal with the possibility that the <& field, which solves the strong-C-P problem, simul- 
taneously accounts for both DE and DM. The angle 6 in eq. (1), as in the PQ model, is turned into 
a dynamical variable, i.e. the phase of a scalar field <£. While 6 is gradually driven to approach 
zero, by our cosmic epoch, in our model <j> gradually increases and approaches m p , yielding DE. 
Residual 6 oscillations, yielding axions, account for DM. 

The main topic of this paper is however the fit of this and other models with WMAP data. 
We compared ACDM, SUGRA dynamical and coupled DE models, as well as a scheme we dubbed 
_1 -model against these data. The last model encloses DM and DE and assumes that they are 
coupled in a non-parametric way, with C{4>) = (j)^ 1 . This last model does not prescribe the origin 
of DE and DM, just as in standard theories of coupled DE. The fit with WMAP constrains A, and 
this constraints agree with the scale range required by the dual-axion model; thus WMAP data 
support the dual-axion scheme. 

The fits of WMAP data to ACDM, uncoupled and constant-coupling SUGRA models, as 
well as to (j)^ 1 SUGRA models, yield similar x 2 's, for all models. At variance from other model 
categories, however, in 4>~ l models CMB data constrain A. This is due to the stronger effects of <fi 
variations on the detailed ISW effect, as they affect both DE pressure and energy density, as well 
as DE-DM coupling. In principle, this strong impact of <j> variation could badly disrupt the fit and 
make (f>~ 1 models significantly farther from data. This does not occur, while the observational A 
range agrees with the dual-axion model at a 2-<r level. 

The success would be complete if the favored range of values of the Hubble parameter (h ~ 0.8- 
1) could be slightly lowered. This range is however obtained just for a SUGRA potential. This 
choice is not compulsory and, moreover, contributions to axion DM due to topological singularities 
were also disregarded. Furthermore, primeval fluctuations were assumed to be strictly adiabatic 
while, in axion models, a contribution from isocurvature modes can be expected. This could 
legitimately affect the apparent position of the first peak in the anisotropy spectrum, so completing 
the success of the model, in a fully self-consistent way. 
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